Excitons with anisotropic effective mass 



On 
G\ 

<d 
in 

(N 

Oh" 

<D 
O 

• i-H . 

CZ5 . 

>V 
^ • 



> ■ 
\o ■ 
m ; 
o . 
o\ ■ 
o ■ 
r» : 

On . 

o : 

>\ 
^ : 

Oh- 



13 



Arno Schindlmayri 

Cavendish Laboratory, University of Cambridge, Madingley Road, Cambridge CB3 OHE, UK 
(Received 25 February 1997, in final form 19 May 1997) 

Abstract. We present a simple analytic scheme for calculating the binding energy of excitons 
in semiconductors that takes full account of the existing anisotropy in the effective mass, as a 
complement to the qualitative treatment in most textbooks. Results obtained for excitons in gallium 
nitride form the basis for a discussion of the accuracy of this approach. 

Zusammenfassung. Wir prasentieren ein einfaches analytisches Verfahren zur Berechnung der 
Bindungsenergie von Exzitonen in Halbleitern, das die vorhandene Anisotropie in der effektiven 
Masse vollstandig miteinbezieht, in Erganzung zu der qualitativen Betrachtung in den meisten 
Lehrbiichern. Ergebnisse fur Exzitonen in Galliumnitrid bilden die Grundlage fur eine Diskussion 
der Genauigkeit dieser Methode. 



I. INTRODUCTION 

The traditional one-particle theory of semiconductors 
as taught in undergraduate courses is that of a mate- 
rial with a finite energy gap separating the highest occu- 
pied from the lowest unoccupied electronic state. In this 
picture the minimum energy for an elementary excita- 
tion is that required to raise a valence band electron into 
the conduction band, and is thus equal to the gap E g . 
However, the electron and the hole, which is created in 
the same process, need not separate completely and can 
instead form a bound pair under the influence of their 
mutual Coulomb attraction M. Such bound electron- 
hole pairs, which transport energy and momentum but 
no charge, are called excitons. They are, in fact, the 
truly lowest elementary excitations of a pure semicon- 
ductor and their influence on the optical properties of a 
material is profound. Most importantly, the occurrence 
of excitons lowers the threshold for photon absorption to 
E g — £"b, where denotes the internal binding energy 
of the electron-hole pair. 

Because of their practical significance excitons feature 
in all textbooks on solid state physics, but the discus- 
sion is usually restricted to some qualitative arguments 
based on formal similarities to the hydrogen atom prob- 
lem, which is modified only by a dielectric constant to 
account for the surrounding medium and by empirical 
effective masses for the electron and hole The latter 
are always assumed to be isotropic, but while there cer- 
tainly are textbook examples such as CdS for which this 
condition is nearly satisfied and hydrogenic absorption 
series have indeed been observed || , in most semiconduc- 
tors the anisotropy in the effective mass is so large that 
it cannot be ignored in a quantitative treatment. In this 
paper we present a variational scheme for calculating the 
binding energy E^ of excitons in realistic materials that 
takes full account of the existing anisotropy. Our priority 
has been to maintain a universally applicable and strictly 
analytic approach suitable for teaching purposes. 



II. VARIATIONAL EXCITON WAVEFUNCTION 

Most semiconductors used for modern electronic de- 
vices crystallize in the diamond (e.g. Si, Ge), zincblende 
(e.g. GaAs) or wurtzite (e.g. GaN) structure, for which 
the constant energy surfaces E(k) about the valence band 
maximum and conduction band minimum arc ellipsoidal 
in shape, yielding distinct longitudinal electron and hole 
effective masses mjj and mj^ along one principal axis, and 

transverse effective masses m c ± and m\ in the plane per- 
pendicular to it. Here we will focus on this geometry, 
although the method is readily generalized. Taking the 
principal axis in the z-direction, the Hamiltonian for the 
relative motion of the electron-hole pair is 



H = 



dx 2 



dy 2 



2/i|| dz 2 



(1) 



where fj,± = m c L m\/(m c L +m^_) and fi\\ = mjjmj^/(m| + 
to!?) denote the reduced transverse and longitudinal effec- 
tive mass, respectively, and e is a suitable dielectric con- 
stant. The anisotropy destroys the spherical symmetry of 
the hydrogen Hamiltonian, yielding a wavefunction with 
a different characteristic localization along the principal 
axis and in the transverse plane. We therefore choose a 
generalization of the hydrogen ground-state wavefunction 
with ellipsoidal symmetry 



ip(x,y,z) = [ — 



1/2 



exp 



(-/V 



x 2 + y 2 + (z/X) 2 (2) 



where the two parameters (3 and A can be varied indepen- 
dently to control the transverse and longitudinal exten- 
sion. The wavefunction, of course, becomes exact in the 
isotropic case = fi± with A = 1 and /3 _1 = aoem/jjL±, 
where a = 0.529 A denotes the Bohr radius and m is the 
free electron mass, but it remains an excellent approxima- 
tion even when the anisotropy is large. A wavefunction 



1 



of the type (£|) was originally proposed in [Q to describe 
shallow donor impurity states in Si and Ge, but treated 
numerically and evaluated only for the mathematically 
distinct case A < 1, reflecting the fact that the longitudi- 
nal effective mass is greater than the transverse one for 
electrons at the bottom of the conduction band in both 
materials. No such restriction will be made in this paper. 

Given the explicit form of the wavefunction (|2|) the 
calculation of the kinetic energy is straightforward, using 
the substitution z' = z/X and a subsequent transforma- 
tion to spherical coordinates. We obtain 



2/xj 



d 2 

ip^-^ipd r 



if 

-P 2 



6/i. 



(3) 



for the contribution in the transverse isotropic plane, and 
similarly 



2/X|| J W dz*- 



6A 2 /Z|| 



(4) 



for the relative motion along the principal axis. To cal- 
culate the potential energy we use the same transforma- 
tion to spherical coordinates. The integrals over the ra- 
dial variable and the azimuth angle are readily evaluated, 
leaving 



1 



-(3 



=-0d 3 r 



sin 9 AO 



87re e J Q ^1 + (A 2 - 1) cos 2 



(5) 



The evaluation of the remaining integral over the polar 
angle depends on the sign of the factor A 2 — 1. If A > 1 
we substitute t = y/X 2 — 1 cos^, otherwise we use the 
substitution t = \f\ — A 2 cos# to obtain an elementary 
integral, which is solved by 



7(A) 



sin 6 AO 



o s/l + (A 2 - 1) cos 2 



arcsinh \/ X 2 — 1 



arcsin yl — 



A 2 



yr^i 2 



for A > 1 



for A < 1. 



(6) 



While we have to make this formal case distinction, we 
emphasize that the energy function is smooth, with both 
branches of I(X) approaching unity in the limit A — > 1. 

Collecting the kinetic and potential contributions, we 
thus obtain the expression 



G 



■ft 2 



2 

M-L 



A M|| 



47rene 



0I(X) (7) 



for the ground state energy, which must be minimized 
with respect to the parameters (3 and A. The kinetic 



term is quadratic in [3 while the potential term is linear, 
so the respective partial derivative is readily performed. 
The condition dE/dj3 = then yields a relation between 
the two parameters at the energy minimum 



ft = 3 



1 



AVii 



/(A) 



(8) 



which, when substituted in (Q), allows us to rewrite the 
energy as a function of A only 



E(X) 



2 y^TTCQeh 



2 

A*_L 



A 2 ^n 



/(A) 2 



(9) 



The energy minimum is found at the stationary point 
for which dE/dX = and through simple mathematical 
rearrangement we can express this condition in the form 



Mi 

M|| 



2A 



1 - XI(X) 



(10) 



which may be solved graphically. The important point 
to note is that the right-hand side of ( |l0| ) is a univer- 
sal function /(A) that does not depend on the material 
properties, so the same plot can be used for all semicon- 
ductors to determine A. In practice, however, the reduced 
transverse and longitudinal effective mass will often not 
differ by more than a factor of three and in this value 
range the function on the right-hand side of ( ]To| ) is accu- 
rately approximated by its lowest-order polynomial term 
/(A) « A 3 . The parameter A is then explicitly given by 



M|| 



1/3 



(11) 



This simplification allows for a very efficient analytic 
calculation of the material-specific binding energy E^, 
which is given by the modulus of the ground state en- 
ergy according to (||). It is still exact in the isotropic 
limit = n±, yielding the correct binding energy 
Eh = RooiTi/([i±e 2 ) where — 13.6 eV is the hydro- 
genic Rydberg energy. The wavefunction itself, required 
for instance to calculate the optical absorption coeffi- 
cient, is given by the original variational expression with 
j3 defined through the relation (§). 



III. NUMERICAL RESULTS FOR GaN 

In order to illustrate the numerical quality of our 
scheme, we now consider the case of gallium nitride as 
an explicit example for excitonic binding energies in a 
realistic semiconductor. In GaN the valence band maxi- 
mum comprises three almost degenerate subbands, giving 
rise to three distinct hole types that can partake in the 
formation of excitons. Conventionally, these are referred 
to as light holes, heavy holes and split-off holes, reflecting 
their different effective masses. Within the same material 
we can study excitons for which the reduced longitudinal 
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TABLE I. Comparison between approximate and numer- 
ically calculated binding energies Eh for excitons in GaN, 
formed with holes of different effective mass m h from the 
three nearly degenerate subbands at the valence band maxi- 
mum. The analytic approximation is accurate even when the 
anisotropy ' s very different from unity. 











E h ( 


meV) 


Hole type 


m\ 


m\ 


A*_l/a*|| 


approx. 


numer. 


Light hole 


0.15 


1.10 


0.483 


15.456 


15.458 


Heavy hole 


1.65 


1.10 


0.959 


24.809 


24.809 


Split-off hole 


1.10 


0.15 


1.805 


18.939 


18.941 



effective mass is greater than the transverse as well as 
excitons for which it is smaller, and test the accuracy of 
the approximate treatment in either case. 

The effective mass of the electron at the bottom of 
the conduction band is mj_ = 0.18 in the transverse and 
mjj = 0.20 in the longitudinal direction, given in units of 
the free electron mass m, while corresponding parameters 
for the three hole types are listed in table |[ All values are 
quoted from ||. Next the calculated anisotropy factors 
M-l/Mh are given. For excitons formed with heavy holes 
this is close to unity, indicating a rather small pertur- 
bation from the isotropic case, but the ratio is substan- 
tially smaller with light holes and larger with split-off 
holes, respectively. The fifth column in table | lists the 
binding energies Eh obtained analytically from (||) with 
A = (/i^/zJn) 1 / 3 . These are compared with results that 
we obtained by exact diagonalization of the Hamiltonian 
(Q) using standard numerical techniques. The applicable 
value of the static dielectric constant is the low-frequency 
limit e = 9.5 of GaN. 

The excellent agreement between the approximate and 
numerical results in all cases confirms the validity of the 
variational wavefunction ( M) a s well as the accuracy of the 
additional simplification ( |ll|) that makes the approach 
strictly analytic. Comparison with experimental data is 
more problematic because it is very difficult to extract 
accurate binding energies from optical spectra. Never- 
theless, the recently published values of 21 meV for exci- 
tons with light and heavy holes and 23 meV for excitons 
with split-off holes in GaN are probably quite reliable 
The discrepancy with the results obtained here is 
due to the underlying model Hamiltonian (|l]), however, 
not the analytic approximations introduced to solve it. 
While the anisotropy in the effective mass is treated ad- 
equately, other important features such as the mixing 
of states at the threefold degenerate valence band max- 
imum or the spatial variation of the dielectric function 
are still neglected. A more involved quantitative scheme 
will also have to incorporate these in order to reproduce 
experimental binding energies for real materials. 
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